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Abstract 

We study the chiral properties of quenched domain wall fermions with several gauge actions. 
We demonstrate that the residual chiral symmetry breaking, which is present for a finite number 
of lattice sites in the fifth dimension (Lg), can be substantially suppressed using improved gauge 
actions. In particular the Symanzik action, the Iwasaki action, and a renormalization group im¬ 
proved gauge action, called doubly blocked Wilson (DBW2), are studied and compared to the 
Wilson action. All improved gauge actions studied show a reduction in the additive residual quark 
mass, TTires- Remarkably, in the DBW2 case m^es is roughly two orders of magnitude smaller than 
the Wilson gauge action at = 2 GeV and Lg = 16. Significant reduction in mj-es is also realized 
at stronger gauge coupling corresponding to a~^ = 1.3 GeV. As our numerical investigation indi¬ 
cates, this reduction is achieved by reducing the number of topological lattice dislocations present 
in the gauge field configurations. We also present detailed results for the quenched light hadron 
spectrum and the pion decay constant using the DBW2 gauge action. 
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I. INTRODUCTION 


Domain wall fermions i I, I, i are expected to provide an implementation of lattice 
fermions with exact chiral symmetry, even at a hnite lattice spacing. To achieve this exact 
symmetry, an inhnite hfth dimension must be introduced and simulations have been done 
to explore the limit of large hfth dimension for both full and quenched QCD I, i 0, S i 
I0| . The hnite size of the hfth dimension, L^, used in numerical simulations, produces a 
small amount of chiral symmetry breaking, which should go to zero in the limit Lg ^ oo. 
In practical implementations the aim is to achieve the smallest chiral symmetry breaking 
possible at a given L^, thus minimizing the cost of the simulation. Further information 
about domain wall fermions and their applications is given in the recent reviews 0 . 0 . 

There have now been several suggestions on how to minimize the computational cost 
of domain wall fermions. An obvious way to achieve this is to make the hve-dimensional 
eigenvectors of the domain wall fermion operator, which for small eigenvalues should be 
localized on the four-dimensional boundaries of the hfth dimension, decay faster in the hfth 
dimension. This reduces the mixing between the opposite chirality modes, which are bound 
to opposite ends of the hfth dimension. Shamir |]^ has calculated the hfth-dimensional 
decay of the eigenfunctions with zero eigenvalues using perturbation theory, suggesting a 
modihcation of the four-dimensional component of the domain wall fermion operator to 
increase the decay. This interesting perturbative result may explain some of the features 
seen in non-perturbative simulations. (Of course, modihcations to the domain wall fermion 
operator may increase the computational cost by more than the reduction in Lg reduces it.) 


Another method of improving domain wall fermions is proposed in |TJ, . The basic idea 
behind these proposals is to project out the zero modes of the four-dimensional Hamiltonian 
describing the propagation in the hfth dimension. As a result, the localization on the 
boundaries of the fermionic light modes is enhanced. 

In this paper we systematically examine a diherent option: the modihcation of the gauge 
action to suppress the hnite Lg explicit chiral symmetry breaking [|^ . Note that in principle 
this is a diherent criteria from improving the gauge action to achieve better scaling, in 
lattice spacing, of physical observables. We will investigate the scaling of observables as 
well, to check that while reducing the explicit chiral symmetry breaking, we do not distort 
the approach to the continuum limit. It is worth noting that methods which improve the 
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domain wall fermion operator, like those suggested by Shamir, and the one investigated here 
are likely independent of each other, so a combination of both techniques may lead to even 
greater efficiency in domain wall fermion simulations. However, as we will see, our approach 
obviates the need of separately treating the near unit eigenvectors of the transfer matrix, 
as gauge conhgurations for which these occur are suppressed. This has also been studied 

m. 


m 


The observation that the gauge action can affect signihcantly the chiral symmetry of 
domain wall fermions is not new. Both the RBC El and CP-PACS [^, ^ collaborations 
have observed that the use of the Iwasaki action 0 substantially improves chiral symmetry 


in quenched simulations. Also in it was observed that the 1-loop Symanzik improved 
gauge action improves chiral symmetry to a lesser degree. Here we extend these results and 
explore the reason behind the observed improvement. 

This paper is organized as follows. In section || we give a brief description of the gauge 
actions under study. In section 0 we introduce the observables used for studying chiral 
symmetry breaking and also present the standard Wilson action results to provide a refer¬ 
ence point. Section 0 contains results for the different actions we studied. We hnd that 
the doubly blocked Wilson (DBW2) action |^, ^ gives residual chiral symmetry breaking 
two orders of magnitude smaller than the Wilson gauge action at comparable lattice spac- 
ings and values of Lg. In section we discuss the dominant mechanism of explicit chiral 
symmetry breaking in domain wall fermions, which we hnd is driven by lattice artifacts, or 
dislocations, at the lattice spacings considered. These dislocations occur as the topological 
charge of the gauge held conhguration changes during Monte Carlo evolution. Given this 
large improvement in residual chiral symmetry breaking and the fact that the DBW2 action 
has not been used before with domain wall fermions, in Section 0 we present results for 
some hadronic observables in order to conhrm consistency with quenched simulations using 
other gauge actions, to check scaling with lattice spacing and to lay a foundation for future 


work 25 


II. PURE GAUGE LATTICE ACTIONS 

As mentioned, we study the chiral properties of quenched domain wall fermions with 
Symanzik, Iwasaki, and DBW2 gauge actions. These actions are built from closed loops of 
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up to six links and provide a sample of typical lattice actions used to improve scaling of 
observables. As a baseline for comparisons we start with the Wilson action which is 
defined by 


Sap] = E 


X,flU 


( 1 ) 




where P\U]x^^u is the real part of the trace of the path ordered product of links around the 
1x1 plaquette in the /i, v plane at point x and f3 ~ l/^fg with qq the bare gauge coupling. 
This is the original non-abelian gauge action introduced by Wilson, which has 0{a‘^) errors 
(a is the lattice spacing). 

To begin, we study the Symanzik one loop improved action where both 0{a'^) and 
0{g‘^a^) errors are removed. This action is defined as 




-f ("^0 + Cl + C2 j 

\ x;fj,<u x;fj,<u<a / 


( 2 ) 


where and C[U]x,fiua denote the real part of the trace of the ordered product of 

SU(3) link matrices along 1x2 rectangles in the p, z/ plane and the p, z/, a, —p, —z/, —a 
paths, respectively. The coefficients cq, ci, and C 2 are computed in tadpole improved one 
loop perturbation theory |^|. For this action and the remaining ones, /3 ~ l/^fg as for the 
Wilson action, but the precise numerical factors differ. 

In addition to the above actions we also studied the Iwasaki action and the DBW2 
action |2^. These actions are both renormalization group (RG) improved actions in a 
truncated, two-parameter space. They can be written down as 

/5 


Sg[U] = 


(1-8 Cl) 


p[u]x,^iu p Cl R[u]x,ij,uj 

X](1<V X^fl^U / 


( 3 ) 


with Cl = —0.331 for the Iwasaki action and ci = —1.4069 for the DBW2 action. In the case 
of the Iwasaki action the coefficient Ci is computed in weak coupling perturbation theory. 
For the DBW2 action ci is computed non-perturbatively using Swendsen’s blocking 


and the Schwinger-Dyson method. QCD-TARO has studied the RG flow in the two 
parameter space of the plaquette and the rectangle couplings and concluded that DBW2 
is a good approximation to the RG flow in this plane at least for a range of coarse lattice 
spacings. 

Although the Iwasaki and DBW2 actions are motivated by the desire to remain on the 
RG trajectory for quenched QGD, the truncation to the explicit form used is an approxi- 
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mation. The accuracy with which these truncated actions preserve the RG trajectory must 
be investigated numerically. Simulations with the Iwasaki action and the DBW2 action 
show improved scaling of the heavy quark potential and the critical temperature for the 
hnite temperature phase transition, compared to the Wilson gauge action. These actions 
serve as useful starting points for studying the effects of the gauge action on residual chiral 
symmetry breaking in domain wall fermions. 


III. EXPLICIT CHIRAL SYMMETRY BREAKING WITH DOMAIN WALL 
FERMIONS 


The central idea behind domain wall fermions is that four-dimensional fermionic states of 
opposite chirality are localized dynamically on opposite boundaries of an extra hfth dimen¬ 
sion. The domain wall fermions are coupled to four-dimensional gauge helds replicated in the 
hfth direction, so the light states can be used to simulate a vector gauge theory like QCD. 
The hve-dimensional fermion action is a generalization of the Wilson fermion action 
with open boundary conditions in the hfth dimension |^. In the free held limit, localization 
of a single fermionic havor on the four-dimensional boundaries occurs if the hve-dimensional 
fermion mass is in the interval (0,2). This interval is shifted when interactions are 
switched on. For an inhnite hfth dimension {Lg —>■ cxd), chiral symmetry of the light states is 
manifest since they have no overlap. Four-dimensional light quark q, q helds are constructed 
from the hve-dimensional fermions T, T by 


q{x) 

= Pl^{x, 0) -h Pr^{x, Lg - 1) 

(4) 

q{x) 

= '^{x,Lg - 1 )Pl + '^{x,0)Pr, 

(5) 


where Pr/l = |(1 i 75 ) are the right-handed and left-handed projection operators. Hence 
a four-dimensional mass term nifqq can be introduced if the hfth dimension boundaries are 
coupled directly with a parameter For hnite Lg explicit chiral symmetry breaking is 

induced by the mixing of the light states which now extend across the hfth dimension. Our 
conventions throughout this paper are the same as those in [§. 

In order to quantify the explicit chiral symmetry breaking induced at hnite Lg, we dehne 
the residual mass {rrires) through the Ward-Takahashi identity |^: 


A,(A“(a;)0(2/)) = 2mf{J^{x)0{y)) + 2{J^^{x)0{y)) + z{6-0{y)), (6) 
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where 


Ls-l 




s=0 




( 7 ) 


is a four-dimensional partially-conserved axial current which is constructed from the hve- 
dimensional conserved vector current, 

. 1 






The flavor matrices are normalized to obey Tr = 5“^, A^f{x) = f{x) — f{x — fi) is a. 
simple hnite difference operator, and the pseudoscalar density J^^x) is 


J^{x) = -T(a:,L,-l)PLt“^(x,0) + T(x,0)PRrT(a:,L,-l). 


(9) 


Note that J^{x) is a four-dimensional pseudoscalar density constructed from helds on the 
boundaries of the hfth dimension. The identity differs from the continuum expression 
by the presence of the J^q{x) term. J^q{x) is analogous to J^{x), but is built from helds in 
the bulk at Ls/2 and Ls/2 — 1. 


J^x) = -T(x, Ls/2 - l)PLp^(a;, L,/2) + T(a;, L,/2 )PrPT(x, L,/2 - 1). (10) 


We refer to this term as the “mid-point” contribution to the divergence of the axial 
current. In the low energy limit the effect of this explicit chiral symmetry breaking is 
described by a simple added residual mass term rures so that in this limit ~ rriresJ^ [H- 
Thus, from the mid-point term we dehne the ratio 


Ex,y('75(l/)^)'75(p 0)) ’ 


( 11 ) 


which, for t greater than some tmin should be independent of t and equal to the residual 
mass, giving 




Ex,y(' 75 g( 2 />^)' 75 (P 0 )) 


T.x,y{J5iy^t)Ji{x, 0 )) 


( 12 ) 


t>tri 


As we will see, in our numerical simulations R{t) is essentially t independent for f ^ 5 and 
the t dependence for t ^ 5 will be discussed in Section To calculate rures, we average 
over a suitable plateau where R(t) is constant. In the subsequent discussion rures serves as 
our basic measure of chiral symmetry breaking. In addition, it is useful to dehne the ratio 


nu]{t) 


'^x,y('^5qiy^ ^)'75 ( 3 ^) 0))[U] 
T.x,y{J^iy^'t)J^{x, 0))[t7] ’ 


(13) 
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and 


(14) 


nu] = 

t 

which are both measures of chiral symmetry breaking on a given gauge conhguration U. 

Before presenting results for the improved gauge actions, we discuss what is known about 
the Wilson action at a~^ ~ 2 GeV {f3 = 6.0). In Fig. ^ we show the residual mass as a 
function of (the data are from [^). While in perturbation theory m^es is expected to 
decay exponentially, as stated in P] the data do not support this. However, its behavior 
can be £t with two exponentials with a rather weak decay in the large Lg limit. Thus, to 
decrease rures by an order of magnitude we need to increase Lg by a large factor, perhaps of 


0 ( 10 ). 

Since rrires is determined by the fifth-dimensional falloff of the boundary states, decreasing 
TTires requires improving the falloff. Analytic arguments have shown that for gauge held 


satisfying a smoothness condition, exponential falloff is assured [p 8 | , . It is expected that 

at weak enough couplings, such a smoothness condition is satished, which is not the case for 
Wilson gauge lattices at /? = 6.0. Since the falloff in the hfth dimension can be related to 
eigenvalues of an appropriately dehned transfer matrix, T, in the hfth dimension, studies 
of the spectrum of the T for Wilson gauge action have been done. They hnd a non-vanishing 
density of unit or near unit eigenvalues of T, showing that undamped propagation in the 
hfth dimension occurs. We will also study the spectrum of T, using gauge conhgurations 
generated with the Wilson, Symanzik, Iwasaki and DBW2 actions. 

The transfer matrix T [^| is dehned by 


l-Ht 

1 + n 


with 


nt = 


2 -|- 




(16) 


being the Hamiltonian for propagation in the hfth dimension and D^{rn) being the four¬ 
dimensional Wilson Dirac operator. Following calculate the eigenvalue spectrum of 

the Hermitian Wilson Dirac operator as a function of M 5 (the so-called spectral 

how). From Eq. one sees that a zero eigenvalue in corresponds directly to 

a unit eigenvalue of the transfer matrix, i.e. the existence of a hve-dimensional mode that 
is not damped in the hfth dimension. In addition, the number of zeros in the spectral how 
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determines the index of the domain wall fermion operator and hence serves as a dehnition 
of topology on the lattice. Thus, if one is working at a hxed value for M 5 and a gauge held 
is generated via Monte Carlo which has a unit eigenvalue of T, an undamped mode in the 
hfth dimension occurs on that conhguration. This conhguration is one where we informally 
say that topology is changing (in the Monte Carlo update). 

When studying the spectral how on a given conhguration, if the how approaches the 
Ms-axis, we expect the left and right domain wall modes to become delocalized leading to 
mixing and attendant chiral symmetry breaking. On the other hand, if there is a large 
vertical gap in the spectral how for values of M 5 we use in our simulations, the chiral modes 
should remain localized on the boundaries. In Fig. ^ we present the spectral how of the 
lowest hfteen eigenvalues for some representative Wilson gauge action conhgurations. Many 
crossings of the M 5 axis are evident and even the modes that do not cross are not far away 
from the axis, compared with the large gap that appears for M 5 < 0.8. Note that M 5 ~ 0.8 
corresponds to the usual critical mass for Wilson fermions where chiral symmetry is restored 
at this gauge couphng(/? = 6 . 0 ). As we will see, this picture leads to a relatively large value 
of rures for the Wilson gauge action, though we emphasize that the chiral symmetry breaking 
is still very small compared to standard Wilson fermions at this gauge coupling. In Fig. ^ the 
ratio r[; 7 ] (t) dehned in Eq. is plotted for the same conhgurations as in Fig. ^ The panels in 
Fig. I and Fig. are in one to one correspondence. In the hgures, rp-^{t) is quite dependent 
on t, with large huctuations occuring over a small range of t. Since we can see multiple 
crossings in the spectral how, which implies undamped modes in the hfth dimension, and 
multiple spikes of rp] {t) it is natural to investigate whether these are diherent manifestations 
of the same phenomena. 

In Fig. ^(a) we present fp] as a function of conhguration number. It is clear that fp] 
huctuates widely, indicating that there are conhgurations with larger chiral symmetry break¬ 
ing and others with relatively small breaking. The number of conhgurations with enhanced 
chiral symmetry breaking is signihcant (~ 50%), consistent with the known result that 
the transfer matrix has an appreciable number of near unit eigenvalues [0. In addition. 
Fig. I suggests a close correlation between conhgurations showing these spikes and those 
with crossings in the spectral how near M 5 = 1 . 8 . 

In order to further examine the nature of chiral symmetry breaking on a given conhgu¬ 
ration we take a closer look at the ratio rp]{t) dehned in Eq. |^. In Fig. ^ and Fig. j^a) 
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we present this ratio for typical Wilson gauge action configurations, again at a~^ ~ 2 GeV. 
As we can see the dominant part of chiral symmetry breaking comes from localized regions 
in time. In particular for the conhguration of Fig. |^(a), the Hermitian Wilson Dirac oper¬ 
ator has two small eigenvalues whose eigenvectors are localized around the peaks of r[[/](f). 
In addition, very localized peaks of a topological charge density, constructed as described 
55[|, are found around the peaks of riu](t). For further discussion see |^. Thus we 


m 


see that the large fluctuations in r[u]{t) as a function of t are due to localized gauge held 
huctuations which are changing the topology of the lattice. This is a crucial observation in 
understanding why improved gauge actions can reduce chiral symmetry breaking for domain 
wall fermions. 

If an improved action can reduce lattice artifact conhgurations which are undergoing 
topology change, then m^es can be reduced. The effect of the gauge action on dislocations 
can be understood by examining its effects on the classical minima of the action i.e. in- 


stantons. Using the results of ||^ we can see that for the Iwasaki and the DBW2 action 
the 0{a?/p^) correction to the action of an isolated lattice instanton is positive, hence in- 
stantons of small size p are suppressed. On the contrary for the Wilson action the 0[a^ jp^) 
correction is negative, consequently the small lattice instantons are enhanced. This suggests 
that for the Iwasaki and the DBW2 actions, gauge conhgurations with very localized concen¬ 
trations of topological charge density are suppressed. If in addition, there is a suppression 
of conhgurations where localized topology change is occuring, there will be a reduction of 
explicit chiral symmetry breaking. In conclusion, conhgurations of non-zero topology do 
not produce large residual chiral symmetry breaking, only conhgurations where topology is 
changing, i.e. where the spectral how has a zero. Suppression of lattice artifact toplogy 
changing conhgurations should decrease rures- 


IV. CHIRAL SYMMETRY WITH IMPROVED GAUGE ACTIONS 

In order to study the ehects of the choice of the gauge action on the residual chiral 
symmetry breaking we performed a series of quenched simulations using the Symanzik, 
Iwasaki and DBW2 actions. In all cases the lattices were 16^ x 32 with inverse lattice spacing 
a~^ ~ 2 GeV. We used the p mass to set the scale but also conhrmed consistency with the 
scale set from the string tension; both yield equal lattice spacings to within a few percent. 
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The mass M 5 was tuned to be optimum with an accuracy of about 5%. Simulations on a 
few conhgurations at several values of M 5 were all that were needed for this determination. 
It turns out that for all actions at a~^ ~ 2 GeV the optimum value is roughly 1.8, except 
for the DBW2 for which it is 1.7. In the free held limit, the optimum value is M 5 = 1[Q. 
The bare quark masses in our study ranged from ruf = 0.010 to 0.060. A summary of the 
simulation parameters is presented in Table |. 

The residual mass was extracted by htting to a constant at large time separations the 


ratio dehned in Eq. |^. Errors are determined by the jackknife method. As it can be seen in 
Fig-i this ratio exhibits a fairly stable plateau at time separations larger than hve or six, 
so we chose a htting range of 7-16 in all cases. All data in this hgure are for Lg = 16 and for 
bare quark mass 0.020. The quark mass dependence of the residual mass is mild as seen in 
Fig.0 Since we have also matched the lattice spacings, it is safe to compare all the actions 
at the same bare quark mass ignoring renormalization ehects. Because the numbers we are 
comparing diher by orders of magnitude these ehects can be safely neglected. In fact the 
multiplicative quark mass renormalization constants have been computed and shown 
to be equal within 5%. In order to eliminate some of the ehects of the remaining small 
mismatch of the lattice spacings, we have plotted the residual mass scaled by the square 
root of the string tension. 

In Fig. ^ we present our measurements of mres for each action for several values of Lg. 
In the case of the Iwasaki action we only performed the measurement at Lg = 16, and our 
result agrees with that of CP-PACS The remaining Iwasaki points are from the CP- 
PACS publication^. As one can see at Lg = 16, the DBW2 residual mass is about two 
orders of magnitude smaller than the residual mass of the Wilson action while the Iwasaki 
residual mass is about an order of magnitude smaller than that of the Wilson action. Finally, 
the residual mass of the Symanzik action is roughly a factor of three smaller than that of the 
Wilson action. In this figure the solid lines represent fits to simple exponentials in all cases 
except the Wilson action where a £t to two exponentials is shown. For the Symanzik data 
a small deviation from the simple exponential £t is observed at Lg = 16 while the Wilson 
action shows a very clear deviation. On the contrary, both the Iwasaki and DBW2 data can 
be fit well with a simple exponential for the same range of Lg. For that reason it is interesting 
to quote a value for the parameter q that Shamir has computed perturbatively 0 - His one 
loop result is that the light fermion wave function x{s) decays exponentially away from the 
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wall, i.e x('S) ~ with q = \. The residual mass also behaves as m^es ~ q^’’■ In the 
case of the Wilson and possibly the Symanzik action, the fact that no good £t to a single 
exponential is obtained may be a signal that m^es scales as a power law [^, and g ~ 1 . 
Such behavior is consistent with the spectral flows observed for the Wilson gauge action. 
For the Iwasaki and DBW2 actions q ~ 0.7 and q ~ 0.6, respectively, which is consistent 
with a gap in the spectral flow at M 5 = 1.7 — 1.8 that is well dehned on most conhgurations. 
We come back to this point in the following where we investigate the spectral flow for each 
gauge action. 

Given the dramatic improvement in rures for the DBW2 action, it is natural to wonder 
whether further improvement is possible. We have explored simulations where the coeffi¬ 
cients of the plaquette and rectangle term in Eq. ^ take on various ratios and found that 
for ratios not far from the DBW2 choice, small lattice spacings could not be achieved. In 
addition, a double peaked distribution of the plaquette values could be found by changing 
the ratio of the plaquette and rectangle coeffecients to be about a factor of 2 different than 
for the DBW2 action. Thus, further dramatic improvement in mres does not seem possible 
with an action which involves only plaquette plus rectangle terms. 


V. TOPOLOGY AND CHIRAL SYMMETRY BREAKING 

In this section we take a closer look at how the different gauge actions affect explicit 
chiral symmetry breaking in domain wall fermions. As mentioned before, in Fig. H(a) the 
quantity fp] defined in Eq. is presented as a function of the conhguration number. The 
large fluctuations (spikes) indicate that there are conhgurations with relatively large chiral 
symmetry breaking and conhgurations with relatively small breaking. The conhgurations 
with large spikes are those for which the transfer matrix in the hfth dimension has a near 
unit eigenvalue, or a corresponding (near) zero eigenvalue of the hermitian Wilson Dirac 
operator. In those cases that we have checked for the Wilson gauge action, a spike is always 
accompanied by a localized (near) zero eigenvector of the the Wilson Dirac operator. In 
addition, the fact that the spectral hows presented in Fig. ^ have so many crossings very 
close to the simulation point M 5 = 1.8 is consistent with the large number of spikes in 
Fig. 1(a). In conhgurations where a spike does not occur, i.e. no crossing close to M 5 = 1.8, 
the chiral symmetry breaking is controlled by the size of the gap of the bulk modes in the 
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spectral flow. Here we are separating the small Wilson Dirac eigenvalues into two groups: 
those that cross zero near M 5 and those that form a more continuum band which we refer 
to as bulk modes. In the case of the Wilson action and conhgurations with no crossings 
close to M 5 , the bulk mode gap is rather small and not very well dehned; thus even on these 
conhgurations the chiral symmetry breaking is relatively large for a given M 5 . 

For the Symanzik action (Fig. ^(b)) the number of spikes is slightly smaller than in the 
case of the Wilson action, and also the number of crossings in the spectral how (Fig. ^ is 
correspondingly reduced. Also, the bulk mode gap is larger. As a result the baseline, or 
level of the troughs between peaks in fp], is lower than in the case of the Wilson action, 
contributing to the reduction in the residual mass. 

The above picture becomes much clearer with the Iwasaki (Fig. §(c)) and DBW2 actions 
(Fig. §(d)). The number of spikes is signihcantly smaller, and the baseline is well dehned 
(especially for the DBW2 action). The typical spectral hows presented in Fig. 0 and Fig. 0 
again support the fact that the Iwasaki action, and to a larger degree the DBW2 action, 
signihcantly suppress the near unity eigenvalues of the domain wall fermion transfer matrix. 
In both cases the gap of the bulk modes in the spectral how becomes signihcantly larger. As 
a result the explicit domain wall fermion chiral symmetry breaking is signihcantly reduced. 

In Fig. 1^ we present the ratio rp]{t) dehned in Eq. |0 for a typical conhguration of each 
action. In all cases it is evident that the dominant contribution to chiral symmetry breaking 
comes from very localized objects, and thus as we argued before, it is not very surprising that 
local 0{a‘^) modihcations of the gauge action can have a very signihcant ehect on explicit 
residual chiral symmetry breaking. 

It is important to recognize that the above mechanism for explicit chiral symmetry break¬ 
ing is related to topology-changing conhgurations (see and references therein for a more 
complete discussion). The connection is made through the index theorem: the domain wall 
fermion operator in the limit > 00 has an index | 0 ] equal to the number of right- 
minus the number of left-handed zero modes, which corresponds to the topological charge 
of the background gauge held conhguration—a quantity which becomes precise in the con¬ 
tinuum limit. This integer depends on the value of M 5 used and is given by the net number 
of crossings in the spectral how of the Wilson Dirac operator as the Wilson mass varies be¬ 
tween a value above the critical Wilson mass and —M 5 . While this index is well-dehned only 
in the limit Lg ^ 00 , our simulations show that the near-zero eigenvectors of the hnite-L^ 
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operator obey the index theorem to a high degree of accuracy |g, In particular, for an 
Iwasaki ~ 2 GeV ensemble, when compared to the topological charge computed using 

I , the index agrees very well. In those cases where 


the smoothing method described in 


the topological charge is not close to an integer, we also hnd a crossing in the spectral flow, 
a spike in r[[/](f), and a complex structure of eigenvectors that is not expected from simple 
chiral symmetry arguments ||38||. If M 5 sits exactly on a crossing, then the index is not de- 
hned, even the limit Lg ^ 00 . A crossing in the spectral flow that occurs away from the 
critical Wilson mass corresponds to a conhguration with indistinct topology. Put differently, 
if the particular gauge held in question is in the midst of changing its topology, which must 
happen if the update algorithm is ergodic and updates the conhguration smoothly, then such 
a gauge held must give rise to a crossing. It is also sensible that such a tunneling from one 
topological sector to another proceeds through local changes in the gauge held which have 
a characteristic size of one to two lattice spacings. In the continuum limit, if the density of 
these dislocations is zero, then all crossings happen at the critical mass and correspond to 
physical topological charge. Thus the index as computed from the spectrum of the domain 
wall operator Dirac operator is well-dehned in this case. 

Consequently, when the Iwasaki action or the DBW2 action is used, the question arises 
whether the topology changes efficiently. We have measured the topological charge using 
the smoothing method described in [P2| , [53[1 ||52|| . Our data are presented in Fig. The 
conhgurations shown in this hgure are separated by 1000 sweeps of Cabibo-Marinari pseudo- 
heatbath with a Kennedy-Pendleton accept/reject step. We can see that there is a 
signihcant slow down in the topological charge fluctuations for the DBW2 action. Both the 
Symanzik and the Iwasaki action also show a mild reduction in the frequency of change of 
the topological charge. Although the problem seems severe for the DBW2 action, we can 
tackle it with brute force. For that reason we have produced a library of DBW2 lattices to 
be used for several domain wall fermion projects. Given the considerable cost of measuring 
domain wall fermionic observables, this higher cost of producing DBW2 lattices at 2 GeV 
is negligible. However, it is clear that this brute force approach will become less practical 
for smaller lattice spacing since topology change is likely to be rapidly suppressed as we 
approach the continuum limit H]. 
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VI. HADRONIC OBSERVABLES FOR THE DBW2 ACTION 


In this section we discuss various hadronic observables calculated with the DBW2 gauge 
action at /? = 0.87 and 1.04 which correspond to a~^ ~ 1.3 and 2 GeV respectively. 

A. The heavy quark potential 


We measure the heavy quark potential as in by fixing to Coulomb gauge and then 
computing the two point correlation function of products of temporal links. More precisely, 


with 


Ae ^W^ = ^(Tr Lt,,T{x)Ll^^^{x + r) ), 

to,X 

to+T 

Lto,T{x) = JJ Ut{x,t), 

t=to 


(17) 


(18) 


and V{r) the heavy quark potential. The potential is extracted by taking ratios of the 
correlation function in Eq. at T and T + 1. The systematics involved in choosing T were 
carefully studied and the optimal T was chosen. For the 1.3 GeV lattices T was 4 while for 
the 2 GeV lattices it was 7. The potential V (r) is fit to 


a 


V{r) = C -h or 


(19) 


The above formula gave very good fits for spatial distances r > a/2. The upper range of r 
was determined by the distance where the error on the potential became unacceptably large. 
The maximum distance used was r = 6 and 7 for the 1.3 and 2 GeV lattices, respectively. 


Figs. [Ij show the heavy quark potential as a function of distance. The results for the 
string tension a and the Sommer parameter ^ vq are tabulated in Table ITTIi These 
results are used in our subsequent discussion of the scaling of hadronic observables. 


B. Simulation and analysis 

For each set of gauge configurations, domain wall fermion propagators are computed with 
two types of sources: 1) a local point source and 2) a Goulomb gauge fixed extended source 
which is either a wall source for (3 = 1.04 or box source with 8^ volume for (3 = 0.87. (We 
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set the source to one at each site inside the box and zero elsewhere.) The local source is 
used for the determination of the decay constants and also the axial current renormalization 
factor Za (/5 = 0.87 only). The extended source is used for all other purposes. 

In Table. we give ruj-es in the chiral limit for the same ensemble of conhgurations used 
for the hadronic observables to be discussed in this section. We have htted mj.es{mf) with 
a linear function of m/ to obtain mres(O) for which the chiral limit of low energy physics is 
dehned as ruf = —mres(O). All data are used to extrapolate to rrif —> 0 for /? = 1.04. On 
the other hand, the largest value rnf = 0.09 for (3 = 0.87 is not used for the extrapolation. 

We take the chiral limit irtf ^ —mres(O) as the physical point for u, d quarks. This 
determines the physical p meson mass rup. With the input rup = 770 MeV, the lattice 
spacing a is determined. The kaon physical point mj, which roughly corresponds to half the 
strange quark mass, for fx and mx* is dehned by mT,{mf = m^)/mp{mf = —rures) = 0.645 
using only degenerate quark masses. We do this procedure for every jackknife sample to 
estimate the error for values at the physical kaon point. 

C. Chiral property of pseudoscalar mass 

Because of the almost exact chiral symmetry of domain wall fermions and the use of the 
quenched approximation, the pion two-point function suffers contamination from topological 
near zero modes, which causes a shift in htted masses from their inhnite volume values[Q. 
The ehect is expected to be inversely proportional to the square root of the volume. Because 
we used diherent physical volumes for the two diherent gauge couplings {V ~ (2.4fm)^ for 
(3 = 0.87, (1.6fm)^ for (3 = 1.04), the size of the ehect on the pseudoscalar mass should be 
diherent in our two lattice ensembles. To study zero mode ehects we examine the pseu¬ 
doscalar mass from two diherent two-point functions. One is the pseudoscalar-pseudoscalar 
correlator (PP) 

(7r“(a;)7r“(0)) = J“(a;) i J“(0)) = -(gr“75g(a;) gr“75g(0)), (20) 

and the other is the correlator of the temporal component the axial-vector current (AA) 

(A“(a;)A“(0)) = gr“757og(0)), (21) 

where q and q are four dimensional quark helds dehned in Eq. |^, As discussed in , the two 
types of correlators suher diherently from topological zero modes. The leading contribution 
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to the pseudoscalar correlator is ~ 1/{rn?j\/V) while it is ~ 1/{rnf\/V) for the axial corre¬ 
lator. Although the relative contribution from the pole compared to the physical one is the 
same for both, the prefactor of the 1/{mf\/V) term is expected to be suppressed Thus, 
the mass extracted from the PP correlator is expected to have a stronger finite volume effect 
from zero modes than the AA correlator. The observed effect on the meson mass calculated 
for light quark masses is to shift it above a linear extrapolation from the region of heavier 
quark mass. 

The pseudoscalar mass extracted from both types of correlators is presented in Table |V]. 
Fig. 1^ shows the pseudoscalar mass squared as a function of m/ for /3 = 0.87. Both values 
of the pion mass are consistent with each other for mj > 0.02. However, at rrif = 0.01 the 
mass extracted from the PP correlator lies above the AA one, outside of their statistical 
errors. Because the axial correlator is expected to have smaller finite volume effects from 
zero modes, we use this correlator for further analysis. 

The linear fit of the pion mass squared in is quite good in the region 0.01 < ruf < 0.06 
as indicated by the x^/dof, which is tabulated in Table ^I[ Note, we are using the same set 
of gauge conhgurations for all values of m/ but employ an uncorrelated fit. One can reliably 
extract the physical kaon mass; however, the fit overshoots the point ruf = —rrires where the 
pion mass should vanish. This is a signal of non-linearity for the pion mass at small mf. 
Instead of a linear function we should employ the quenched chiral log formula with the 
constraint that the pion mass vanishes a.t rrif = —rrij-es, 


ml = a^{mf + m^es) 


(a^{mf + m,es) 
'^log -- 

V ^^QxPT 


+ b{mf + mres)^ 


( 22 ) 


where we have introduced the quadratic term in addition to the expression used in As 
is discussed in Ref. Eq. has a form suggested by chiral perturbation theory. Note 
that the explicit chiral perturbation theory formula, e.g. Eq. 90 of Ref. [^, involves three 
parameters a, S and Kq^pt and additional terms. However, within our parameter range 
and numerical accuracy these additional terms have the effect of adding an undetermined, 
mj-independent constant to the expression within the square brackets. Our choice of zero 
for this unknown constant represents a rescaling of the parameters 5 and Aq^pp from those 
that appear in the chiral perturbation theory prediction. 

The data fit this formula well with a reasonable value of the chiral log coefficient 6 as listed 
in Table [VII[ , where the scale is set as Aq^pp = 1 GeV. If the coefficient of the quadratic 
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term b is not set to zero in the fit, a somewhat larger value of 6 results. This is because 
the two terms tend to cancel each other. Note that the value of b is the same order as a^r. 
Ultimately, a proper covariant £t with reliable should be used to distinguish the two fits. 
At present we do not have enough statistics to do this. 

Fig. p!B| shows the pion mass squared for (3 = 1.04. The results for the larger value of mj 
agree with each other, but deviations appear as m/ decreases. In fact, pronounced upward 
curvature of the data, which was not observed for the axial correlator for (3 = 0.87, result in 
large values of x^/dof listed in Table |V| for both correlators. The lines in the Fig. |l^ show 
the results of the fit excluding the lightest point ruf = 0.01, which makes x^/dof smaller. 
Thus we can assume that the pion mass for ruf > 0.015 does not suffer much from zero 
mode effects. As this £t reproduces the data in the range 0.015 < ruf < 0.04 very well, we 
use it for the interpolation of the kaon point. Also using this region oi ruf, one can £t data 
with the quenched chiral log without the quadratic term for x^/dof = 0.5(1). The resulting 
6 is consistent to that for (3 = 0.87. Since our range of for (3 = 1.04 is not wide enough 
to disentangle the quadratic term and the log term, we do not list a result for non-zero b. 

D. Hadron spectrum 

We list the results for the vector meson and nucleon masses in Table 0. Figs. and 
p!8| show the vector meson and nucleon masses as functions of m/ for (3 = 0.87 and 1.04 
respectively. Physical nucleon, rho, and K* masses are indicated on the figure. 

Fig. |T^ presents the rho meson mass versus lattice spacing squared, both normalized 
with tq. We also plot the results obtained with the Wilson gauge action from ref. 0. We 
have selected values obtained at (/?,!/, L^) = (6,16^, 24), (5.85,12^, 20), (5.7, 8^, 48). These 
lattices have almost the same physical volume, and Ls is the largest available in each case. 
We observe consistency between DBW2 and Wilson actions. The flatness of the data reflects 
the small size of the scaling violation. 

Fig. gg is a scaling plot for the nucleon and K* masses normalized by the p meson mass. 
Again good scaling is observed. The K* appears to be lighter than experiment, which is 
consistent with other quenched simulations. For the nucleon mass, there is an observed 
discrepancy between the results with Wilson type fermions and staggered fermions (see 
comparison by S. Aoki ^4|). The former gives a lighter nucleon mass while the latter is 
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consistent with the experiment. The nncleon mass for domain-wall fermion is slightly larger 
than experiment for the lattice spacings we examined. We need to perform simnlations for 
larger physical volnme at a~^ ~ 2 GeV, as well as simnlations at smaller lattice spacings 
to do the continnnm extrapolation needed to compare with conventional fermions. For 
the comparison within the domain-wall fermions, given the statistics and the fact that the 
physical volnmes of onr ensembles are not the same, we can not say if the DBW2 action 
exhibits better scaling than the Wilson gange action thongh it seems that the scaling is at 
least as good. 


The J parameter, which is introdnced in ref. to examine the effect of qnenching on 
the mass spectrnm, is dehned by 

rim 

(23) 


drrip 

J = rriK*- 


dm^ ’ 

where mK*{mf) is evalnated at rnp{rnf)/mn{mf) = mK*/mK = 1-8 by dehnition. The slope 
of rUp with respect to is determined with the data in 0.01 < ruf < 0.06 for jS = 0.87 
and 0.015 < rrif < 0.04 for f3 = 1.04, nsing pion mass determined with AA. An approximate 
phenomenological valne calcnlated from the experimental mass spectrnm is. 


mK*-mp 

J = rriK *—^- TT = 0.48. 


m 


K 


rat 


(24) 


Onr lattice resnlts for the J parameter, listed in Table |V111| and plotted in £g. |^, show a 
signihcant difference from the phenomenological valne. However, these resnlts are consistent 
with other qnenched results (see a summary given by T. Kaneko |^). Note that we have 
used degenerate quark masses which could explain some of the difference between our result 
and experiment, but the largest source of the discrepancy is likely to be due to quenching. 


E. Pseudoscalar decay constants 

The pseudoscalar decay constants are calculated from the amplitude of the point-point 
two-point function of the temporal component of local axial vector current, 

= (25) 

^ X 

for f S> 1. This simple local current is not partially conserved, unlike Ap (Eq. |^) which 
obeys the Ward-Takahashi identity (Eq. so the local current receives a multiplicative 
renormalization. We can calculate the renormalization factor from the ratio of the 
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correlation functions of the two currents. Here we employ a ratio designed to remove some 
of the lattice spacing error H. 


W) = - 


C{t + 1/2)+ C{t-1/2) ^ 2C{t+l/2) 


(26) 


2L{t) L{t)+L{t + l)_ 

where C and L are the correlators of the pseudoscalar density with the partially conserved 
and local axial currents respectively, 


C(« + l/2) = ^{.4S(J,«)P“(0)), 

X 


(27) 

(28) 


The hrst and second terms on the r.h.s. of Eq. PB| remove the 0(a) scaling error, and suppress 
the 0{a‘^) error in the sum. The value of is determined by htting TZ{t) to a constant for 
the interval Q < t < 26, where TZ{t) is flat to a very good approximation. The results are 
given in Table ^ and plotted in Fig. An estimate of Za with only the hrst term is also 


plotted to demonstrate that the complex ratio in Eq. actually works. Although linear 
dependence on my in TZ remains for both cases, it is very small for the full expression of 
Eq. We use the value Za obtained in the chiral limit my —> —m^es for the calculation of 


the pseudoscalar decay constant, which is given in Table 

Another method to calculate the decay constant is to use the Ward-Takahashi identity 
to derive its relation with the pseudoscalar to vacuum matrix element of the pseudoscalar 
density: 


A 


m; 


= (0|J5“|7r,p = 0). 


(29) 


(my+mres) 2 

This matrix element is determined from the PP correlator. Results with each correlator for 
= 0.87 are shown in Fig. ^ For both methods we use the results of m^ extracted from 
the AA correlator with the extended source quark propagator since it is our most precise 
determination. 

Consistency between both results shows that we have good control over the chiral sym¬ 
metry breaking through our determination of mres- Also, possible zero mode effects, which 
influence the two correlators differently, appear to be small at least for my > 0.02, which 
is consistent with the absence of zero mode effects in the pion mass over the same range of 
my. 


Fig. |2J is the same plot but for /3 = 1.04. The results from the two determinations 
are also in good agreement. Results for decay constants at the physical points are given 
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in Table and plotted in Fig. ^ against lattice spacing squared. They are consistent 
with those reported for the Wilson action p. We observe good scaling both for and 
fx- The result for is in good agreement with the experimental value. However, the ratio 
fx/f-K appears to be smaller than the experimental value 1.22, which is an expected effect 
of quenching 


VII. CONCLUSIONS 


In this paper we have demonstrated that the DBW2 action signihcantly improves the 
chiral properties of domain wall fermions. The main reason for this improvement is the 
suppression of gauge conhgurations which support unit eigenvalues of the transfer matrix 
in the hfth dimension and hence allow signihcant mixing of the light chiral modes that are 
localized on opposite boundaries of the fifth dimension. These problematic configurations 
are also those which occur as the topological charge of the gauge field changes during Monte 
Carlo evolution. One key to improving the domain wall fermion chiral symmetry is to use 
an improved gauge action which suppresses these small dislocations. This suppression works 
so well in the case of the DBW2 action that at a~^ = 2 GeV and Lg = 16 the residual chiral 
symmetry breaking is roughly two orders of magnitude smaller compared to the Wilson 
action case and therefore is completely negligible. Even at strong coupling (a“^ 1.3 GeV) 

rures is about three times smaller in physical units than for the Wilson action at 2 GeV. In 
both cases the value of Lg is 16. 

Besides the suppression of these small topological dislocations associated with zero cross¬ 
ings of the spectral flow of the four-dimensional Wilson Dirac operator, we have also observed 
an increased gap in the spectral flow. Gonsequently, the light boundary mode wavefunctions 
decay faster in the fifth dimension. For the DBW2 action this leads to a residual mass 
dependence on Lg proportional to 0.6^®. This dependence is close to Shamir’s perturbative 
prediction |]^. Approaches based on the proposals made by Shamir in may be effective 
in further reducing this perturbative baseline. Work along these lines is currently underway. 

In the second part of this paper we presented results for some quenched hadronic observ¬ 
ables obtained with the DBW2 gauge action. Our conclusion is that these observables scale 
very well with a, i. e. the good scaling of domain wall fermions seen in quenched simulations 
with the Wilson gauge action is preserved [Q. 
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Using these improved actions, we also observed that the topological charge of the gauge 
helds evolves much more slowly using standard Monte Carlo algorithms. In future simula¬ 
tions with smaller lattice spacings, improved algorithms will be needed to efficiently sample 
topology. We note, however, that this is a generic feature of all lattice calculations which is 
not specific to the DBW2 action. 
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TABLE I: Simulation parameters for each gauge action tested. The p mass, mp, is given in the 
chiral limit for the largest Lg in each case. As usual (3 ~ l/g^, where go is the bare gauge coupling, 
M 5 is the five-dimensional fermion mass, Lg the size of the 5th dimension, m/ the bare quark mass, 
and a is the string tension computed from the heavy quark potential. 


Action 

/3 

Ms 

Lg 

ruf 

rup 


Wilson p 

6.00 

1.8 

12-24 

0.015-0.040 

0.404(8) 

0.227(6) ID 

Symanzik 

8.40 

1.8 

8-16 

0.020-0.060 

0.411(14) 

0.2278(18) II 

Iwasaki 

2.60 

1.8 

16 

0.020-0.060 

0.415(13) 

0.231(6) m 

DBW2 

1.04 

1.7 

8-16 

0.020-0.060 

0.399(11) 

0.2246(16) 

TABLE 11: The residual 

mass rures at a ^ Rs 

2 GeV for the actions tested. In the construction 

of this table, : 

for the Symanzik action we used 51 configurations, for the Iwasaki 45, and for the 

DBW2 89. 







ruf 

Lg 


Symanzik 

Iwasaki 

DBW2 

0.020 

8 


3.04(5)10-3 


- 

7.54(5)10-'^ 

0.020 

12 


8.2(4)10-^ 


- 

9.92(20)10-3 

0.020 

16 


3.3(3)10-“^ 

1.4(4)10-^ 

1.60(5)10-3 

0.040 

8 


2.90(4)10-3 


- 

7.49(5)10-^ 

0.040 

12 


7.4(3)10-“^ 


- 

9.9(4)10-3 

0.040 

16 


2.73(24)10-^ 

1.2(4)10-'^ 

1.56(3)10-3 

0.060 

8 


2.82(3)10-3 


- 

7.50(8)10-^ 

0.060 

12 


6.95(23)10-^ 


- 

1 . 00 ( 6 ) 10 -^ 

0.060 

16 


2.44(18)10-^ 

1.15(27)10-^ 

1.565(23)10-3 


24 














TABLE III: Parameters and resulting scales for the DBW2 gauge action used for the spectrum 
analysis. The quantity Op is the lattice spacing determined by the p meson mass taken from 
Table. |VII1. The jackknife errors are quoted for the string tension cr and the Sommer parameter 


ro- 


/3 

Ms 

Ls 

statistics 


^■o 

a-i(GeV) 

0.87 

1.8 

16 

100 

0.324(6) 

3.58(4) 

1.31(4) 

1.04 

1.7 

16 

405 

0.2246(16) 

5.24(3) 

1.98(3) 


TABLE IV: The residual mass for the DBW2 gauge action calculated on the same configurations 
used to evaluate the hadronic observables. The values for ruf = 0 have been obtained from a linear 
fit in nif. The value for ruf = 0.09 at /3 = 0.87 is excluded from the fit. 



(3 = 0.87 


P = 1.04 

mf 

Tflxes 

ruf 

TTZres 

0.01 

5.44 (23) 10“^ 

0.01 

1.80 (9) 10-5 

0.02 

5.16 (22) 10“^ 

0.015 

1.80 (11) 10-5 

0.03 

4.84 (21) 10“^ 

0.02 

1.77 (11) 10-5 

0.04 

4.55 (19) 10“^ 

0.025 

1.74 (11) 10-5 

0.05 

4.30 (18) 10“^ 

0.03 

1.71 (10) 10-5 

0.06 

4.08 (16) 10-^ 

0.035 

1.69 (8) 10-5 

0.09 

3.52 (13) 10-^ 

0.04 

1.67 (7) 10-5 

0 

5.69 (26) 10“^ 

0 

1.86 (11) 10-5 


25 
















TABLE V: Hadron masses computed using the DBW2 gauge action. The superscripts PP and 
AA refer to the pseudoscalar and axial vector correlation functions, respectively. All masses are 


obtained with degenerate quarks. 

/3 

ruf 

PP 


rrip 

niN 


0.01 

0.2248 (25) 

0.2179 (31) 

0.607 (22) 

0.790 (75) 


0.02 

0.2997 (19) 

0.2966 (23) 

0.640 (17) 

0.871 (23) 


0.03 

0.3603 (16) 

0.3590 (20) 

0.662 (11) 

0.921 (14) 

0.87 

0.04 

0.4128 (15) 

0.4118 (19) 

0.685 (8) 

0.975 (10) 


0.05 

0.4601 (14) 

0.4589 (18) 

0.709 (6) 

1.021 (8) 


0.06 

0.5037 (13) 

0.5021 (17) 

0.732 (5) 

1.067 (7) 


0.09 

0.6192 (13) 

0.6170 (15) 

0.803 (4) 

1.197 (6) 


0.01 

0.1794 (22) 

0.1759 (21) 

0.413 (6) 

0.546 (13) 


0.015 

0.2098 (17) 

0.2075 (18) 

0.424 (4) 

0.575 (9) 


0.02 

0.2377 (15) 

0.2359 (16) 

0.435 (4) 

0.602 (7) 

1.04 

0.025 

0.2631 (13) 

0.2617 (15) 

0.447 (3) 

0.628 (6) 


0.03 

0.2868 (12) 

0.2857 (14) 

0.4586 (29) 

0.652 (5) 


0.035 

0.3090 (12) 

0.3081 (13) 

0.4705 (26) 

0.676 (4) 


0.04 

0.3300 (11) 

0.3293 (12) 

0.4825 (24) 

0.699 (4) 


TABLE VI: Results from fitting to the linear function 

fitting range. 

Co + ciuif. The column m/ 

shows the 

/3 

correlator 

nif 

Co 

Cl 


dof 

0.87 

PP 

0.01-0.06 

0.0090 (14) 

4.057 (29) 

2.8 (1.2) 

4 

0.87 

AA 

0.01-0.06 

0.0063 (16) 

4.090 (37) 

0.17 (35) 

4 

1.04 

PP 

0.01-0.04 

0.0056 (9) 

2.566 (21) 

3.4 (8) 

5 

1.04 

PP 

0.015-0.04 

0.0047 (9) 

2.597 (19) 

1.10 (25) 

4 

1.04 

AA 

0.01-0.04 

0.0044 (9) 

2.585 (21) 

2.6 (7) 

5 

1.04 

AA 

0.015-0.04 

0.0035 (9) 

2.615 (20) 

0.86 (23) 

4 
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TABLE VII: Fit results for the quenched chiral log contribution to the pion mass determined from 
the AA correlator using eq. ^ The symbol indicates that the parameter is constrained to be 
zero. 


/3 

ruf 

Q-tt 

b 

6 


dof 

0.87 

0.01-0.06 

4.04 (5) 

- 

0.031 (14) 

1.6 (1.5) 

4 

0.87 

0.01-0.09 

3.40 (23) 

6.4(1.6) 

0.107 (38) 

0.19 (25) 

4 

1.04 

0.015-0.04 

2.583 (28) 

- 

0.049 (14) 

2.0 (5) 

4 


TABLE VIIL Results from fitting the vector meson mass to a linear function, cq + ciuif. Values 
at the physical points and the J parameter are also listed. 


/9 

mf 

Co 

Cl 


dof 

rup 

mK*/mp 

J 

0.87 

0.01-0.06 

0.590 (19) 

2.37 (28) 

0.14 (50) 

4 

0.589 (19) 

1.138 (11) 

0.377 (37) 

1.04 

0.01-0.04 

0.389 (6) 

2.34 (12) 

0.08 (20) 

5 

0.388 (6) 

1.136 (5) 

0.387 (16) 


TABLE IX: Results from fitting the baryon mass to a linear function, cq + ciruf. Values at the 
physical point are also listed. 


/3 

ruf 

Co 

Cl 

X^ 

dof 

ruN 

0.87 

0.01-0.06 

0.783 (27) 

4.76 (44) 

0.5 (1.6) 

4 

0.780 (27) 

1.04 

0.01-0.04 

0.502 (12) 

4.96 (26) 

0.5 (8) 

5 

0.502 (12) 
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TABLE X: The axial current renormalization factor and pseudoscalar decay constants. 


/3 

mf 

Za 

fAA 

JPS 

fPP 

J PS 


0.02 

0.78199 (37) 

0.1084 (31) 

0.1078 (34) 


0.03 

0.78404 (31) 

0.1121 (27) 

0.1110 (27) 

0.87 

0.04 

0.78612 (28) 

0.1165 (26) 

0.1145 (24) 


0.05 

0.78824 (26) 

0.1208 (26) 

0.1192 (24) 


0.06 

0.79042 (25) 

0.1247 (26) 

0.1237 (24) 


0.01 

0.84142 (17) 

0.0703 (22) 

0.0697 (21) 


0.015 

0.84191 (14) 

0.0716 (17) 

0.0720 (17) 


0.02 

0.84244 (12) 

0.0732 (14) 

0.0743 (15) 

1.04 

0.025 

0.84300 (11) 

0.0751 (12) 

0.0767 (14) 


0.03 

0.84358 (10) 

0.0771 (11) 

0.0790 (13) 


0.035 

0.84417 (9) 

0.0790 (11) 

0.0813 (12) 


0.04 

0.84478 (9) 

0.0809 (10) 

0.0835 (12) 


TABLE XL The axial current renormalization factor at ruf = 

—nires and the pseudoscalar decay 

constants at physical points in [MeV]. 






/3 

ZA{mf = -rUres) 

fAA 

J 77 

fPP 

J 77 

fAA 

JK 

fPP 

JK 

KP/ff-' 

fP/fP 

0.87 

0.77759(45) 

130.4(6.7) 

129.0(7.3) 

148.9(5.2) 

147.3(5.4) 

1.142(26) 

1.141(30) 

1.04 

0.84018(18) 

130.8(4.9) 

129.0(5.0) 

147.4(3.3) 

149.7(3.6) 

1.139(24) 

1.118(25) 














0.020 



10 20 30 40 50 


FIG. 1: The residual mass normalized by the square root of the string tension (a) as function of 
the Lg for the Wilson gauge action at /3 = 6.0 (results are from |^). The solid line is a fit to a 
double exponential. At large Lg the decay is rather weak. 
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FIG. 2: Spectral flows of the Hermitian Wilson Dirac operator 'y^Dw for four typical gauge con¬ 
figurations generated with the Wilson gauge action at (3 = 6.0. There are many crossings in the 
neighborhood of M 5 = 1.8 which induce explicit chiral symmetry breaking for domain wall fermions. 
The size of the would-be gap in the region of the five-dimensional fermion mass M 5 ~ 1.8 is also 
relatively small compared to the obvious gap above the critical Wilson mass M 5 ~ 0.8. Both effects 
enhance mixing of the light domain wall fermion modes and hence the value of the rrires- 
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r[u](t) r[u](t) 




FIG. 3: The ratio of Eq. |I^ for four Wilson gauge configurations. The bare quark mass is 0.02 
and M 5 = 1.8. The configurations used are the same as those in Fig. and the panels are in one 
to one correspondence with the panels in Fig. || 
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FIG. 4: The quantity defined in Eq. 14 vs. configuration number for (a) the Wilson gauge 
action at /? = 6.0, (b) the Symanzik gauge action at /3 = 8.4, (c) the Iwasaki gauge action at 
(3 = 2.6 , (d) and the DBW2 gauge action at (5 = 1.04. All four cases correspond to a~^ ss 2 GeV. 
The large spikes seen in the Wilson case are significantly reduced for the Iwasaki action and almost 
eliminated for the DBW2 action. These spikes corresponded to zero eigenvectors of and are 

a significant source of chiral symmetry breaking for domain wall fermions. 
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FIG. 5; The ratio of Eq. |T^ for a single configuration of (a) the Wilson gauge action at /3 = 6.0, 
(b) the Symanzik gauge action at /3 = 8.4, (c) the Iwasaki gauge action sX (5 = 2.6, (d) and the 
DBW2 gauge action at /3 = 1.04. The spikes are quite localized in Euclidean time t. Examination 
of the eigenvectors of the (domain wall fermion or Wilson) Dirac operator confirms the zero modes 
are localized in space as well. 


33 








s 



t 


FIG. 6: The ratio defined in Eq. 11 at a 


-1 


2 GeV. The fancy squares correspond to the Wilson 


gauge action, the diamonds to Symanzik, the squares to Iwasaki, and the octagons to DBW2. The 


bare quark mass in all cases is 0.020 and Lg = 16. 
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FIG. 7: The residual mass at a ^ ~ 2 GeV as a function of the bare quark mass. The octagons 
correspond to DBW2, the squares to Iwasaki, and the diamonds to Symanzik. In each case Lg = 16. 
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FIG. 8: Dependence of the residual mass on the size of the 5th dimension at ~ 2 GeV. The 
octagons correspond to DBW2, the squares (GP-PACS [^) and diamond (RBG Q) to Iwasaki, 
the bursts to Symanzik, and the fancy squares to Wilson. All but the Wilson action fit a simple 
exponential decay reasonably well. Note the Iwasaki results use different gauge field ensembles 
at each value of Lg. In the case of the Wilson action, the results are fit to a double exponential 
function. 
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eigenvalue eigenvalue 



FIG. 9: The same as Fig. but for the Symanzik gauge action at /3 = 8.4. The number of crossings 
in the neighborhood of M 5 = 1.8 appears slightly smaller and the would-be gap slightly larger than 
in the Wilson case. 
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FIG. 10: The same as Fig. but for the Iwasaki gauge action at /9 = 2.6. The number of crossings 
neighborhood of M 5 = 1.8 is significantly smaller and the gap clearly larger than in the Wilson 
case. 
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eigenvalue eigenvalue 




FIG. 11: The same as Fig. but for the DBW2 gauge action at (3 = 1.04. There are no crossings 
neighborhood of M 5 = 1.8 and the gap is quite large, roughly comparable to the gap at the 
corresponding mass above the critical Wilson mass. Note that this is also true for the region 
beyond the next critical Wilson mass, M 5 ~ 2.3, where the four flavor Wilson fermion doublers 
become light. 
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conf # 


FIG. 12: The configuration history of the topological charge Q. The configurations are separated 
by 1000 sweeps of Cabibo-Marinari pseudo-heatbath with a Kennedy-Pendleton accept/reject step. 
At the top is the Symanzik action with lattice size 16^ x 32 and coupling (3 = 8.4. In the middle 
is the Iwasaki action with size 16^ and /3 = 2.6. At the bottom is the DBW2 action with size 
16^ X 32, and /3 = 1.04. All cases correspond to roughly the same scale, a~^ ~ 2 GeV. The DWB2 
action, which suppresses configurations with small instantons, shows a significant reduction of the 
tunneling between topological charge sectors. 
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FIG. 13: The heavy quark potential for DBW2, f3 


0.87. The solid line denotes the fit to Eq. IS 


from which the string tension is determined. 
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FIG. 14: The same as Fig. 13 but for (3 = 1.04 
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FIG. 15: The quenched pion mass squared as a function of for (5 = 0.87. Masses extracted from 
the pseudoscalar (PP) and axial vector (AA) correlation functions are shown. They agree quite 
well, except for the point at rrif = 0.01 where there is a small difference outside of statistical errors. 
The physical volume is roughly (2.4 fm)^, so contamination from topological zero modes which can 
induce such splittings (see text) is suppressed. The line is a fit to a simple linear function. The 
extrapolation slightly overshoots the expected chiral limit point rrif = —nires (see inset), suggesting 
a quenched chiral log. 
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FIG. 16: The same as Fig. |]^ but for j3 = 1.04. 
zero mode effects in the masses are visible, and 
limit point ruf = —rures- 


Here the physical volume is roughly (1.6 fm)^, so 
the simple linear fit clearly overshoots the chiral 
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FIG. 18: Vector meson and baryon masses as functions of mj for (3 = 1.04. The line is a linear fit 
to the data. The physical K* mass is also shown. 



FIG. 19: Scaling of the p meson mass with the lattice spacing set by the Sommer parameter tq. 
The Wilson data are from ref. 
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FIG. 20: Scaling of the nucleon and K* masses with lattice spacing. The experimental mass of 
nucleon and K* are shown on the vertical axis. The Wilson data are from ref. |^]. 



FIG. 21: Scaling of the quenched J parameter with lattice spacing. The experimental point is 
shown on the vertical axis. The large discrepancy with experiment is due to quenching. 
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FIG. 22: The local axial current renormalization factor as a function of mj for /3 = 0.87. The lines 
correspond to fits to Eq. (lower) and the first term only (upper). The difference arises from 
lattice spacing errors that are proportional to ruf. 
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FIG. 23: The quenched pseudoscalar decay constant as a function of for /3 = 0.87. The results 
from the two techniques described in the text agree quite well. The data are fit to simple linear 
functions, and extrapolated (pion) or interpolated (kaon) to the physical points (filled symbols). 
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FIG. 24: The same as Fig. but for (3 = 1.04. Although the results agree within statistical 
errors, the agreement is not as good as that at /3 = 0.87. Since the physical volume is significantly 
smaller in this case, we expect the difference is due to the presence of topological zero modes. 
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FIG. 25: A scaling plot for the pion and kaon decay constants. The DBW2 results appear to 
preserve the good scaling that was observed with the use of the Wilson gauge action (data from 
Ref. i) 
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